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We review some recent results in digital geometry obtained by using a combinatorics on words ap- 
proach to discrete geometry. Motivated on the one hand by the well-known theory of Sturmian words 
which model conveniently discrete lines in the plane, and on the other hand by the development of 
digital geometry, this study reveals strong links between the two fields. Discrete figures are iden- 
tified with polyominoes encoded by words. The combinatorial tools lead to elegant descriptions of 
geometrical features and efficient algorithms. Among these, radix-trees are useful for efficiently de- 
tecting path intersection, Lyndon and Christoffel words appear as the main tools for describing digital 
convexity; equations on words allow to better understand tilings by translations. 

1 Introduction 

The expansion of computers has led in the last few decades to several breakthrough in technological 
achievements. Among these, digital imaging is increasing its spread and is extensively used in a wide 
range of applications such as image synthesis, remote sensing, medical image processing to cite a few. 
Developed mostly by the engineering world, digital geometry has led to the discovery (sometimes a 
rediscovery) of new results about discrete sets, concurrently to the design of new algorithms tools, and 
enriched the broader field of discrete geometry. 

Combinatorics on words has imposed itself as a powerful tool for the study of discrete, linear, and 
non-commutative objects that appear in almost any branches of mathematics, and discrete geometry is not 
an exception. Traditionally, digital geometry works on characterization and recognition of discrete ob- 
jects using an arithmetic approach or computational geometry. However combinatorics on words provide 
some useful tools and efficient algorithms for handling discrete objects. Lothaire's books |[20l |2T1 l22l 
constitute the reference for presenting a unified view on combinatorics on words and many of its appli- 
cations. 

As mentioned by Klette and Rozenfeld in their survey on digital straightness lITSl 

"Related work even earlier on the theory of words, specifically, on mechanical or Sturmian 
words, remained unnoticed in the pattern recognition community" 

there was a need for new investigations where combinatorics on words would enrich the classical Eu- 
clidean approach of digital geometry. 

We revisit some classical problems in discrete geometry from this new point of view. For our purpose 
the discrete plane is identified with the square grid Z x Z. 
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2 Preliminaries 

We refer to Lothaire |[20l for the basic terminology and notation about words on a finite alphabet £/. 
It includes the empty word e, length, , conjugate, factor, prefix, suffix, proper factor, free monoid si* , 
morphism, antimorphism, occurrences, palindrome, period, power, primitive, reversal. The set of all 
factors of w is denoted by Fact(H'), those of length n is Fact„(H') = Fact(H') fl Pref(w) is the set of 
all prefixes of w, and the set of its palindromic factors is Pal(L). If w = pu, with |w| = « and \p\ =k, 
then p^^w = w[k..n — 1] = m is the word obtained by erasing from w its prefix p. The class of a word 
w is denoted [w]. Every word contains palindromes, the letters and e being necessarily part of them. 
This justifies the introduction of the function LPS : jz/* — > s^* which associates to any word w its longest 
palindromic suffix LPS(w). Given a total order < on si, the lexicographic ordering is defined as usual. 

Lyndon words Introduced as standard lexicographic sequences by Lyndon in 1954, Lyndon words 
have several characterizations (see flOl). We shall define them as words being strictly smaller than any 
of their circular permutations. 

Definition 1. A Lyndon word I € si^ is a word such that I = uv with u,v (z s/^ implies that I < vu. 

Note that Lyndon words are always primitive. The most important result about Lyndon words is the 
following unique factorization theorem (see Lothaire |[20l Theorem 5.1.1). 

Theorem 1. Any wordw ^ admits a unique factorization as a sequence of decreasing Lyndon words: 



where ni > 1 and is a Lyndon word, for all i such that \ < i <k. 

There exist several algorithms for factorizing a word w = wiW2 • • ■ w„ into Lyndon words and the 
more efficient are linear . An elegant one was invented by Duval ||161 . It works by reading from left to 
right, with at most 2n comparisons of letters (see also [31 1, Section 7.4). Another one, uses the concept 
of suffix standardization of the word w, and builds a suffix array of w, which may be computed in linear 
time llT4l . Then the Lyndon factorization of w is obtained by cutting w just before each left-to-right 
minimum of its suffix-array. 

A quadtree witli a radix tree structure for points in tlie integer plane (||9l[T0l) Let B = {0, 1} be 

the base for writing integers. Words in B* are conveniently represented in the radix order by a complete 
binary tree (see for instance |[T9l l22|), where the level k contains all the binary words of length k, and the 
order is given by the breadth-first traversal of the tree. To distinguish a natural number x G N from its 
representation we write jc G B*. The edges are defined inductively by the rewriting rule x ^ x -O+x - l, 
with the convention that and 1 are the labels of, respectively, the left and right edges of the node having 
value X. This representation is extended to B* x B* as follows. As usual, the concatenation is extended 
to the cartesian product of words by setting for G B* x B*, and (a,j8) G B x B 



/I 



'ni,n2 lilt 
1 '2 



h>h> ■■■>lk 



(1) 



(jr,j)-(a,j8) = {x-a,y- 



Let X and y be two binary words having same length. Then the rule 



{x,y) {x-Q,y- 



Q) + {x-0,y\) + {x-\j- 



0) + {x-\,y\) 



(2) 



defines a 5^' 



{N,R), sub-graph of ^ = {N ,R,T 



), such that : 
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(i) the root is labeled (0,0); 

(ii) each node (except the root) has four sons; 
(ill) if a node is labeled {x,y) then \x\ = |j|; 

(iv) edges are undirected, e.g. may be followed in both directions. 

By convention, edges leading to the sons have labels from the ordered set {(0,0), (0, 1), (1,0), (1, 1)}. 
These labels equip the quadtree with a radix tree structure for Equation ^ implies that {x' ,y') is a son 
of {x,y), if and only if 

(x ,y) = {2x+a,2y + ^), 

for some (a, j8) G B x B. Observe that any pair {x,y) of nonnegative integers is represented exactly once 
in this tree. Indeed, if = |j| (by filling with zeros at the left of the shortest one), the sequence of pairs 
of digits (the two digits in first place, the two digits in second place, and so on) gives the unique path in 
the tree leading to this pair. Of course the root may have up to three sons since no edge labeled (0, 0) 
starts from the root. 




(0,0) (0,1) (1,0) (1,1) / (0,0) (0,1) (1,0) (1,1) 




Figure 1: The point (2, 1) with its neighbors. 



Neighboring links ll9j[T0l Given {x,y) G 1?, a point {x' ,y') is an e-neighbor of {x,y) if there exists 
e such that {x',y') = {x,y) + e = {x + ei,y + e2)- 

We superpose on G' the neighboring relation given by the edges of T (dashed lines). More precisely, 
for each elementary translation e G each node (?) = {x,y) is linked to its £-neighbor (2) + £, when it 
exists. If a level k is fixed, it is easy to construct the graph 

^#« = (A^«,/?«,r«) 

such that 

(i) if (xj) G A^W, then \x\ = \y\ = t, 

(ii) the functions A^t^^' ■— > N x N ^ B* x B* are injective; 

(iii) /?(^' is the radix-tree representation : (B<'^ x B<'=) x (B x B) A B-'^ x B-*^; 

(iv) the neighboring relation is r(*=) C x (B x B) x A^. 
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Note that the labeUng in Fig. [T]is superfluous: each node represents indeed an integer unambiguously 
determined by the path from the root using edges in R; similarly for the ordered edges. Moreover, if a 
given subset M C N x N has to be represented, then one may trim the unnecessary nodes so that the 
corresponding graph '^m is not necessarily complete. 

Recall that adding 1 to an integer jc S B*^ is easily performed by a sequential function. Indeed, every 
positive integer can be written x = uVO^, where / > 1, 7 > 0, with u G {e} U jlB*^^'^^^^ -0} . In other 
words, 1^ is the last run of I's. The piece of code for adding 1 to an integer written in base 2 is 



If j y^O then Return ul'OJ-^l; 

else If M = £ then Return 1 • 0' ; 

else Return m • 0^^ • 1 • 0'; 

end if 

end if 



where 0^' means to erase a 0. Clearly, the computation time of this algorithm is proportional to the 
length of the last run of I's. Much better is achieved with the radix tree structure, where, given a node 
(z), its father is denoted /((?)), and we write f{x,y) or f{x,y) if its label is {x,y). The following technical 
lemma is a direct adaptation to B* x B* of the addition above. 

Lemma 2. Let G^'^'^ be the complete graph representing M-^ x M-^for some k>\,e&.^, and (z) = {x, y) 
be a node ofN'^. If one of the four conditions holds: 

(i) £ = and jc[k] = 0, (ii) £ = 2 and jc[k] = 1, 
(iii) £ = 1 and j[k]=0, (iv) £ = 3 and j[k] = 1, 

then f{@) =f{® + £). Otherwise, /(@) + £ = /(@ + £). 

The process is illustrated for case (i) in the diagram ^ ••• \ 

on the right where the nodes (101 10, •) and (101 11, •) „ 

1011,0101 [ ^ ^ ^ ^ iioo.oioT] 

share the same father while fathers of neighboring / 
nodes 



(0,1) (1.1) (0,1) 



(•,01011) and (..OlOll) 

are distinct but share the same neighboring relation. 



10110,01011}-"^ 10111,01011^'"^ 11000,01011 



A representation for paths in the square grid Here, we encode paths with the so-called Freeman 
chain code\23\ based on the alphabet ^ = {0,1,2,3}, considered as the additive group of integers 
mod 4. Basic transformations on ^ are rotations p' : x + i and reflections Oi : x ^ i — x, which 
extend uniquely to morphisms (w.r.t concatenation) on ^* . Given a nonempty word w G ^* , the first 
differences word A(w) € ^* of w is 

A(w) = (W2-Wi)-(W3-W2)---(W„-W„_i). (3) 

One may verify that if z S then A{wz) = A(h')A(w„zi)A(z). Words in ^* are interpreted as paths 
in the square grid, so that we indistinctly talk of any word w G ^* as the path w. Moreover, the word 
w := p^(w) is homologous to w, i.e., in direction opposite to that of w (Figure |2ll. A word u G J^* may 
contain factors in = {02,20,13,31}, corresponding to cancelling steps on a path. Nevertheless, each 
word w can be reduced in a unique way to a word w', by sequentially applying the rewriting rules in 
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Figure 2: (a) w = 01012223211. (b) A(w) = 1311001330. (c) w = 33010003232. 



{m 1-^ £ I M G "^j. The reduced word w' of w is nothing but a word in = ^* \ ^*'^ . The turning 
number^ of w is defined by 3^{w) = (|A(w')|i - |A(w')|3) /4. 

A path w is closed if it satisfies |w|o = \w\2 and = |w|3, and it is simple if no proper factor of w 
is closed. A boundary word is a simple and closed path, and a polyomino is a subset of I? contained in 
some boundary word. It is convenient to represent each closed path w by its conjugacy class [w], also 
called circular word. An adjustment is necessary to the function ^ , for we take into account the closing 
turn. The first differences also noted A is defined on any closed path w by setting 

A([w]) = A(w) • (wi -w„), 

which is also a closed word. By applying the same rewriting rules, a circular word [w] is circularly- 
reduced to a unique word [ w'] . If w is a closed path, then the turning number^ of w is 

^w) = [w]) = (|A( [w'DIi - |A( [w'DIa) /4. 

It corresponds to its total curvature divided by 271. Clearly, the turning number ^{ [w] ) of a closed path 
w belongs to 2, (see El HI). 

The convex hull of a finite set of points The lexicographic order < on points of or I? is such 
that {x,y) < {x' ,y') when either x < x' or x = x' and y < y'. The convex hull of a finite set S of points 
in is 

the intersection of all convex sets containing these points and is denoted by Conv(5). S being 
finite, it is clearly a polygon in the plane whose vertices are elements of S. The upper convex hull of S, 
denoted by Conv+(5), is the clockwise oriented sequence of consecutive edges of Conv(5) starting from 
the lowest vertex and ending on the highest vertex. The lower convex hull of S, denoted by Conv^(5'), 
is the clockwise oriented sequence of consecutive edges of Conv(S) starting from the highest vertex and 
ending on the lowest vertex. 

3 The Daurat-Nivat relation [ffi ^ 

We recall from Daurat and Nivat ifTSl that a discrete set £ is a subset of and an element (/, j) € E 
corresponds to a unit square with vertices {i± ± One sets P1/2 = (5) 2) and the salient and 
reentrant points are defined as follows: 

Definition 2. (Daurat and Nivat |[T5l ) Let E be a discrete set. Then 

(i) A comer is a couple (M,N) when M S /'1/2 and N G and M — fi is in 

(ii) A comer (M,N) is salient ifN G E and M is the extremity of two edges of the border ofE which 
are also sides of the square N + [— ^, ^]^. 

'in (3l|6l, the authors introduced the notion of winding number of w which is AS/'{w). 
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(iii) A comer (M,N) is reentrant j/N G £ and M is the extremities of two consecutive edges of the 
border E which are not sides of the square N + [— 2, 5]^- 

Definition 3. (Daurat and Nivat[i5|) The multiset of border-salient (resp. border-reentrant) points of a 
discrete set E, denoted Sb{E) (resp. Rb{E)} is the multiset whose support is included in P\j2 and such 
that for any M S P1/2. number mult5g(£)(M) {resp. multsg(£)(M)j is the number ofN such that 
(M,N) is a salient comer (resp. reentrant comer). 

In other words, in their terminology modulo a translation by (5,5), a point M on the boundary of 
a polyomino P is salient (see Figure |3] (a)) if it belongs to the intersection of two consecutive sides of 
a square belonging to P. The point M is reentrant (see Figure [3] (b)) if it is the intersection of two 
consecutive edges of the contour of P which are sides of a square not belonging to P. 




(a) (b) 
Figure 3: Salient and reentrant points in tlie terminology of Daurat and Nivat 

In our framework, a salient point of the boundary w of a polyomino corresponds to a left turn (a 1 in 
A(w)) and a reentrant one to a right one (a 3 in A(w)), provided the traversal is done in a counterclockwise 
manner. The Daurat-Nivat |[T5l states that the S salient and R reentrant points in every polyomino are 
related by the formula 

S-R = 4. (4) 



N 




Figure 4: A region and its four extremal points. 



The four extremal points are defined by the coordinates: W is the lowest intersection with the left 
side of the bounding rectangle Q, N the leftmost intersection with the top side, E the highest intersection 
with the right side, and S the rightmost intersection with the bottom side. Note that the four extra salient 
points can be canonically identified as W, S, E and N. Then, the Daurat-Nivat relation |[T5]| may be 
restated as 

Proposition 3. The turning number of a boundary word w is ^w) = ±1. 

Note that a boundary word w is positively oriented (counterclockwise) iff its turning number is 
0(w) = 1. 
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Remarks 1. This rather elementary geometrical property is fundamental in proving a result about 
deciding whether a polyomino tiles the plane by translation or not. 

2. This result can be extended easily to take into account hexagonal grids in which case the alphabet 
must be extended to 6 letters, in which case the Daurat-Nivat relation becomes (see S |5l for more 

details) 

S-R = 6. (5) 

3. The statement above includes closed paths not necessarily simple. It cannot therefore be used to 
determine whether a closed path is simple or not. 



4 Path Intersection m W 

Many problems in discrete geometry involve the analysis of the contour of discrete sets and many prob- 
lems are solved by using linear algorithms in the length of the contour word. However, most of the time 
it is assumed that the path encoded by this word does not intersect itself. Checking non intersection 
amounts to check if a grid point is visited twice. Of course, one might easily provide an iff{n\ogn) al- 
gorithm where sorting is involved, or use hash tables providing a linear time algorithm on average but 
not in worst case. The underlying principle of the algorithm is to build a graph W = {N,R, T) where N 
is a set of nodes associated to points of the plane, R and T are two distincts sets of oriented edges. The 
edges in R give a quadtree structure on the nodes while the edges in T are links from each node to its 
neighbors. 



The Algorithm First, we assume that the path is coded by a word w starting at the origin (0,0), and 
stays in the first quadrant N x N. This means that the coordinates of all points are nonnegative. Note that 
in N X N, each point has exactly four neighbors with the exception of the origin (0,0) which admits only 
two neighbors, namely (0, 1) and (1,0), and the points on the half lines (x,0) and (0,^) with x,y >\ 
which admit only three neighbors. 

Now, assume that the node (x,j) exists and that its neighbor {x+ \,y + 0) does not. If |jif| = |j| = fc, 
then the translation (x,y) + (1,0) is obtained in three steps by the following rules: 

1. take the edge in /? to /(jc,j) = {x[\..k-l\,y[l..k-l\); 

2. take (or create) the edge in T from f{x,y) to (?) = f{x,y) + (1,0); 

3. take (or create) the edge in R from (z) to (2) • (0,j[^]). 

By LemmalU we have ®-{0,y[k]) = (x+ 1,3^+0), so that it remains to add the neighboring link {x^y) 
{x+ \,y + 0). Then, a nonempty word w G is sequentially processed to build the graph and we 
illustrate the algorithm on the input word w = 0011. 

• Initialization: one starts with the graph consisting of the 
node (0,0) marked as visited. For convenience, the non- 
visited nodes (0,1), (1,0), and the links from (0,0) to its 
neighbors are also added. This is justified by the fact that 
the algorithm applies to nonempty words. 
Since (0,0) is an ancestor of all nodes, this ensures that ev- 
ery node has an ancestor linked with its neighbors. The cur- 
rent node is set to (0,0) and this graph is called the initial 
graph ^g. 



0,1 



ool 




(0,1) 



(1,0); 



1,0 



Edges in R 

Edges in T 

I I Visited 
I I Non-visited 
Current node 



Figure 5: Initial graph Sfg. 
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• Read w I = 0: this corresponds to the translation (0, 0) + (1 , 0). A neighboring 
link labeled starting from (0,0) and leading to the node (1,0) does exist, so 
the only thing to do is to follow this link and mark the node (1,0) as visited. 
The current node is now set to (1,0), and this new graph is called 



• Read W2 = 0: this time, there is no edge in % labeled starting from (1,0). 
Using the translation rules above, we perform: 

(1) go back to the father /( 1 , 0) = (0, 0) ; 

(2) follow the link to (1,0); 

(3) add node (2,0) ~ (1,0) • (0,0) = (10,00). 

Then an edge from (1,0) to (2,0) with label is added to T. Finally the node 
(2,0) is marked as visited, and becomes the current node. 

• Read wj, = 1: this amounts to perform the translation (2,0) + (0, 1). Since the 
edge to /(2,0) is labeled by (0,0), we know that the second coordinate of the 
current node (2,0) is even. Therefore, (2, 1) and (2,0) must be siblings, that is 
/((2,0) + (0, 1)) = /((2,0)). What we need to do then is : 

(1) go back to the father /(2,0) = (1,0); 

(2) follow the edge 1 if it exists; 

Since it is does not exist, it must be created to reach the node (2, 1) ~ (10,01) = 
(1,0) • (0, 1). Again an edge from (2,0) to (2, 1) with label 1 is added, (2, 1) is 
marked as visited and is now the current node. 



• Read W4 = 1: since /((2, 1)) has no neighboring link labeled by 
1, recursion is used to find (or build if necessary) the node corre- 
sponding to its translation by 1. This leads to the creation of the 
node (1, 1) ~ (0,0) • (1, 1) marked as non-visited. Then, the node 
(2,2) ~ (1, 1) • (0,0) is added, marked as visited, and becomes the 
current node. Note that a neighboring links between (1,0) and (1,1), 
(2, 1) and (2,2) are added in order to possibly avoid searches. 



0,0 



(ai) /;(i,o); 



0,1 1 1,0 1 

Figure 6: Graph ^^o- 



0,0 



(ai)^/;'(i,o)'; 



0,1 



1,0 



(0,0) 



Figure 7: Graph ' 



0,0 

(0^1) /;'(i,o)'; 

1,0 



?00- 



0,1 



/ (0,0) 

iZ 



(0,1) 

2,0h-^H2j1 
Figure 8: Graph Wooi- 



0,0 



((XI) 
6, 



(1,0) 



(1,1) 



0,1 





1,0 


u 




/ \ 




/ (0,0) {0,1} 


1 / 




\ 


2,0 


--^^ 2,1 - 



(0,0) 



Figure 9: Graph %oii- 



This algorithm is linear and we refer to fTOl for details about the complexity analysis which is rather 
involved. 



5 Digital convexityllUlTl 

The notion of convexity does not translate trivially, and detecting if a discrete region of the plane is 
convex requires a deeper analysis. There are several (more or less) equivalent definitions of digital 
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convexity, depending on whether or not one asks the digital set to be connected. We say that a word 
w is digitally convex if it is the boundary word of a finite 4-connected subset S of I? such that S = 
Conv(X)nZ2. Given such a region: 



W is the lowest on the Left side; 
N is the leftmost on the Top side; 
E is the highest on the Right side; 
S is the rightmost on the Bottom side; 
So that w = W1W2W3W4. 



We say that a word wi in ^* is NW-convex iff there are no integer points between the upper convex hull 
of the points {<^WiU)} j=\...\wi\ and the path w,-. Then we have 

Wi is convex <;=^ p'^'(w,) is NW-convex. 

Clearly, the convexity of w requires the convexity of each w,- for / = 1,2,3,4, and we have the following 
obvious property. Observe that if for some /, w,- contains more than 2 letters, that is if Alph(p'^^(H',)) ^ 
{0, 1}, then w is not digitally convex. 

Theorem 4 {fT\). A word v is NW-convex iff its unique Lyndon factorization I'^H'^ ■ ■ ■ is such that all 
li are primitive Christoffel words. 

For example, the Lyndon factorization of the word v 
1011010100010 is 

v = (1)1 -(011)1 •(01)2.(0001)^ -(0)1, 

where 0, Oil, 01, 0001 and are all Christoffel words. 




This result leads to a fast optimal algorithm for checking digital convexity of a boundary word. It is 
based on the linear time algorithms for computing the Lyndon factorization of the contour word and for 
the recognition of Christoffel factors which digital line segments. By avoiding arithmetical computations 
the algorithm is much simpler to implement and much faster in practice (see |7| for more details). 

It is worth noting that many results about Sturmian words have been obtained by using geometrical 
properties. This close relation between the two domains raises a number of combinatorial problems 
such as enumeration of convex words of given length [301. Another interesting property is the factorial 
closure of digitally convex words: while a geometrical proof is rather easy to obtain, it begs for a purely 
combinatorial proof. A formal language characterization is also a challenge. 



6 Tilings 

Tilings appeared as one of the archetypes of the close relationship between art and mathematics, and are 
present in human history under various representations. The beautiful book of Griinbaum and Shephard 
ifTTl contains a systematic study of tilings, presenting a number of challenging problems (see also Q 
for related work). For instance, the problem of designing an efficient algorithm for deciding whether a 
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given polygon tiles the plane becomes more tractable when restricted to polyominoes, that is, subsets of 
the square lattice whose boundary is a non-crossing closed path. Indeed, while a sufficient condition 
is provided by the Conway criterion in L32 |, one of its consequences is that the only objects that tile 
the plane by translation in two directions are generalizations of parallelograms and parallel hexagons, 
hexagons whose opposite sides are equal and parallel (see |32| p. 225 for more details). One may 
consider these tiles are continuous deformations of either the unit square or the regular hexagon. 



Tiling the plane by translating a single polyomino ( llT3l [H) Such a polyomino is called exact in 
||33]| and Beauquier and Nivat [1] characterized them by showing that their boundary word satisfies the 
equation b{P) = X Y Z X Y Z, where at most one of the variables is empty and where W = p^(W). 
This condition is referred to as the BN-factorization. An exact polyomino is said to be a hexagon if 
none of the variables X, Y , Z is empty and a square if one of them is so. While decidability was 
already established in |[33]| . recently, a i?{n\og^ n) was designed for deciding if a word w € ^ tiles the 
plane by translation HllTTl. It uses several data structures that include radix-trees, for checking that 
w is a closed non crossing path [10], and suffix-trees for building a tricky algorithm for checking the 
BN-factorization [,13,1. 



Square tiles In this case there is a linear algorithm for deciding if a word is a square fSl. Moreover, 
it turns out that a square has at most two distinct square tilings ll26l l27l . which means that the BN- 
factorizations are distinct and therefore not conjugate. Motivated by the attempt of characterizing the 
square tiles, the study of equations of the form ABAB = XYXY lead to a subset of solutions in bijec- 
tion with the BN-factorizations. It turns out that the solutions are strongly related to local periodicity 
involving palindromes and conjugate words |[25l 129^ 1. They also give some insight on their shape since 
palindromes represent symmetric sides. An interesting fact is that some infinite families of such tiles, 
namely the Fibonacci tiles have a fractal characteristic and are connected with some problems in number 
theory |[28i . while Christoffel words yield another infinite family of squares ll24ll . 



Remark Again a number of problems are raised by these connections about the enumeration and gen- 
eration of tiling polyominoes. It turns out that the prototype of a square is the basic one cell polyomino, 
and it suffices to replace each pair of sides by homologous paths. It remains only to check that theses 
paths do not intersect, which is achieved in linear time thanks to the optimal algorithm we have. For the 
generation of double square tiles, the prototype is the cross polyomino, and their generation is a bit more 
involved. The case of hexagon grids deserves also some attention HI. 




S. Brlek 
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